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p ; Abstract 

• An aperiodic (low frequency) spectrum may originate from the error term in the mean value of an arithmetical function such as Mobius 

m : function or Mangoldt function, which are coding sequences for prime numbers. In the discrete Fourier transform the analyzing wave is 

periodic and not well suited to represent the low frequency regime. In place we introduce a new signal processing tool based on the 
Ramanujan sums c q (n), well adapted to the analysis of arithmetical sequences with many resonances p/q. The sums are quasi-periodic 
versus the time n of the resonance and aperiodic versus the order q of the resonance. New results arise from the use of this Ramanujan- 
£Nj ' Fourier transform (RFT) in the context of arithmetical and experimental signals. 
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0. INTRODUCTION 
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"In this age of computers, it is very natural to replace the continuous with the finite. One 
thinks nothing about replacing the real line R with a finite circle (i.e., a finite ring Z/qZ) and 
similarly one replaces the real Fourier transform with the fast Fourier transform" [1]. 

In this paper our claim is that the discrete Fourier transform (and thus the fast Fourier 
transform or FFT) is well suited to the analysis of periodic or quasi periodic sequences, but 
fails to discover the constructive features of aperiodic sequences, such as low frequency 
noise. This claim is not new and led to alternative time series analysis methods such as 
Poincare maps [2] (i.e. one dimensional return maps of the form x n+l = f(x n ) or more 
general multidimensional maps), fractal or wavelet analysis methods [3] and autoregressive 
moving average (ARMA) models [4] to mention a few. These methods appeared in 
diverse contexts: turbulence, financial, ecological, physiological and astrophysical data. For 
stochastic sequences such as l/f electronic noise only small progress was obtained thanks to 
these techniques [2]. 

Here we introduce still another approach by considering the time series as an arithmetical 
sequence, that is a discrete sequence x(n), n= 1 • • -t, in which generic arithmetical functions 
(such as o(n) the sum of divisors of n, <j)(q) the number of irreducible fractions of 
denominator q, the Mobius function jLt(n), or the Mangoldt function A(n) may be hidden. 

Recently we published a number of papers which emphasize the connection between 
frequency measurements and arithmetic [2], [5], [6]. The standard heterodyne method, which 
compares one oscillator of frequency f(n) at time n to a reference oscillator of frequency 
f , leads to irreducible fraction p i jq i of index i given from continued fraction expansions of 
v = f(n) / f and beat signals of frequencies F(n) = fnq^V—pJq^. Jumps between fractions 
of index z, i ± 1, z ± 2 • • ■ were clearly identified as a source of white or l/f frequency noise 
in such frequency counting measurements [5]. A phase locked loop was characterized as 
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well, leading to a possible relationship between 1/ f noise close to baseband and arithmetical 
sequences of prime number theory [6]. 

We introduce Ramanujan sums as a new signal processing tool for these experimental 
files. In contrast to the discrete Fourier transform in which the basis functions are all roots 
of unity (Sect. 1), the Ramanujan-Fourier transform (RFT) is defined from powers over the 
primitive roots of unity (Sect. 2). We provide a table of known RFT's and emphasize the 
newly discovered connection between 1/ f noise and arithmetic. In this context considered in 
Sect. 3, a modified Mangoldt function plays a central role, since it connects the Golden 
ratio found in the slope of the FFT to the Mobius function found in the structure of the 
RFT. Concerning the experimental files considered in Sect. 4, galactic nuclei are promising 
candidates for the application of the new method. 



1. THE DISCRETE FOURIER TRANSFORM 



The discrete Fourier transform (DFT) or its fast analogue (the FFT) is a well known signal 
processing tool. It extends the conventional Fourier analysis to sequences with finite period q 
(for the FFT one takes q = 2 1 , with / integer). 

In the DFT one starts with the roots of unity of the form exp(2?7T-), p — 1 . . .q and the 

signal analysis is performed thanks to the n th power 

p 

e„{n) = exp(2?7T— n). (1) 

q 

(In the mathematical language one says that e p (n) is a character of G = Z/qZ, it is a group 
homomorphism from the additive group G into the multiplicative group of complex numbers 
of norm 1). The DFT of the time series x(n) is defined as 

q 

x(p) = (2) 

n=\ 

and they are a number of relations such as the inversion formula 

1 q 

= -J\x(p)e P (n), (3) 

q , 

i p=l 

the Parseval formula (conservation of energy) 

fWn)| 2 =^X>(p)| 2 , (4) 

77=1 Q P =l 

the orthogonality relations between the characters 

n=l < 

and the convolution formula 

x*y = xy, (6) 

where * means the convolution. Generic discrete Fourier transforms are given in table 1. 
In particular, as it is well known, an oscillating signal e,(n) = exp(2i7t^-n) of frequency II q 
transforms to a line at p = I in the DFT spectrum. Inversely a line at n = I in the time series 
transforms to an oscillating signal exp(—2in^p) of frequency l/q. 
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Table 1. 



x(n) 


x{p) 


1 




e,(n) 


q b){p) 


8,(n) 






cos(2np/q) 


Lq(-n) 


L q (n)L q {-\) 



A Gaussian transforms to a Gaussian through the Fourier integral. Not so well known 
is that the role of the Gaussian is played by the Legendre symbol in the context of the DFT 

[1],[7]. Let us define the Legendre symbol L q (n) = 0^) for an odd prime q as follows 

if q divides n, 
+ 1 if n is a square modulo q 

(x 2 = n (mod q) has a solution) , 
— 1 otherwise. 



(7) 



There are a number of relations attached to the Legendre symbol 

A, 



n 2 (mod q) . 



= (-l) V 



q\fP\ p-l q-l 

— ) [ — \ — (~ l)" 5 " for distinct odd primes p,q, 

pj \qj 



SI, 

The invariant relation for the DFT on Z/qZ is 

L q (-n) = g L q (n) with g = L q (-1). (8) 
The Fourier coefficient at position n equals the coefficient of the original sequence up to a 
constant factor g = LJ-l) = ( £ ) exp(2in E ) and g 2 = (-l^q. 

2. THE RAMANUJAN FOURIER TRANSFORM 

Ramanujan sums c q (n) are defined as the sums of the n th powers of the q th primitive roots of 
the unity [8], [9] 



c q (n) 



£ exp(2^n), (9) 
p=i q 

(P,9)=l 

where (p,q) = 1 means that p and q are coprimes. It may be observed that the c q (n) are the 
sums over the primitive characters e p (n). The sums were introduced by Ramanujan to play 
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x q = —f-rA v (x(n)c q (n)). (13) 



the role of basis functions over which typical arithmetical functions x(n) may be projected 

GO 

x ( n ) = £-*<7 c <?( n )- (1°) 

q=\ 

It should be observed that the infinite expansion with q — > °° reminds the Fourier series 
analysis, rather than the discrete Fourier transform that is taken with a finite q. As a 
typical example the function o(n) (the sum of divisors of n) expands with a RFT coefficient 

o q = ^^, that is 

, s 7T 2 n ri (-l) n 2cos(2/zti;/3) 2cos(/ztt/2) 

a ( n )=-6-{ l+ V~ + — — 

For functions x(n) having a mean value 

A v (x) = ]imWx(n), (12) 

one obtains the inversion formula 
1 

<KqY 

More general formula have also been derived [10]. In the rest of the paper the coefficient x q 
given in (13) will be we called the Ramanujan-Fourier transform or RFT. It follows from a 
number of important relations. There is the multiplicative property of Ramanujan sums 

C qq'( H ) = c q( n ) c q>( n ) if te'tf') = X - ( 14 ) 

and the orthogonality property 

qq' 

n=l 
1 

c q{ n ) = otherwise, (15) 

«=l 

which reminds us (5). It is relatively easy to evaluate Ramanujan sums from basic functions 
of number theory. Let us denote (q,n) the greatest common divisor of q and n. Using the 
unique prime number decomposition of q and n 

q = ]Jqf i (q t prime), 

n = ri n f* ("^ prime) , (16) 

k 

(17) 

one gets the number <j)(q) of irreducible fractions of denominator q, also called Euler totient 
function 

^(^)=^ri( 1 --)' 

q\ 

and a coding of prime numbers from the Mobius function \x{n) which is defined as 

if n contains a square fi k > 1 , 
lifn=l, 

(-l)*ifn is the product 1 j 

of k distinct primes. 



Ramanujan sums are evaluated from[l 1] 

^ n) = J')*!S). (20) 

Note that for (q,n) = 1, c q (n) = n(q). The first values are given from 



c x = 1; c 2 = -1,1; c 3 = -1,-1,2 

c 4 = 0, -2,0,2; c 5 = -1,-1,-1,4--- (21) 

where the bar indicates the period. For instance c 3 (l) = — 1, c 3 (2) = — 1, c 3 (3) =2, c 3 (4) = 
— 1... Some generic Ramanujan-Fourier transforms are given in table 2. In the table the 



Table 2. 



x(n) 


Xq 


a(n) 
n 


X 2 1 

6 q2 


m 

n 


6 fl{q) 
X 2 <l> (q) 


b{n) - ^ n)A{n) 




C{n) 





function 2 (?) generalizes Euler function 

^2^)=^ 2 n( 1 -^)- 

In the Mangoldt function A(n) is defined as 

According to Hardy and Littlewood (1922) the number of pair of primes of the form p,p + h 



is 



7C h (x)~C(h) 



X 



ln z {x) 



with 



C(h) 



2C 2^-p\h p-2 
0, 



! — if h odd 
if h even. 



(24) 



(25) 



where p > 2 is a prime, and the notation p|/z means p divides h. The parameter C 2 ~ 0.660... 
is the twin prime constant. It was recently conjectured [ ] that this problem of prime pairs is 
also related to an autocorrelation function from the Wiener-Khintchine formula 



A v (b(n)b(n+h)) = C(h). 



(26) 



3. LOW FREQUENCY NOISE FROM ARITHMETICAL FUNCTIONS 



The idea which subtends our new signal processing is that experimental signals may hide 
arithmetical features. It is thus very important to master the low frequency effects due to 
generic arithmetical functions such as Mobius function, Mangoldt function and so on [5], [6]. 
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Figure 1. The normalized summatory Mobius functionM(f)/? 1 / 2 . 



100 




1 10 100 1000 

Fourier Frequency f 

Figure 2. The power spectral density (FFT) of the normalized summatory Mobius function 
M(t)/t x l 2 in comparison to the power law 1 / f 2 (dotted line). 




resonance q 



Figure 3. The Ramanujan-Fourier transform (RFT) of the normalized summatory Mobius 
function shown in Fig.l. 

3.1. On the summatory Mobius function 
Let us consider the summatory function 

MO) = £ n(n) = 0(M £ ), whatever e. (27) 

n=l 

The asymptotic dependance assumes the Riemann hypothesis [ ]. The normalized summatory 
function M(t)/t 1 / 2 is shown in Fig. 1. The corresponding power spectral density is in Fig. 2; 
it looks like the FFT of a random walk since the slope is close to —2. 

The RFT of M(t)/t x l 2 is shown in Fig. 3. There is no known formula for it, but it shows 
a signature with well defined peaks which is reminiscent of the function fi(q) /<j>(q) shown 
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below in Fig. 8. 

3.2. Results related to the Mangoldt function 

Riemann hypothesis can also be studied thanks to the summatory Mangoldt function 

x i f{t)=Y d A{t) = t{\ + e w {t)). (28) 

The error term represented in Fig. 4 can be expressed analytically from the singularities (the 
pole and the zeros) of the Riemann zeta function [5]. Fig. 5 shows the FFT of the error term 
£y(t): it roughly behaves as 1// noise. 




Figure 4. Error term in the Mangoldt function A(«). 




1 10 100 1000 

Fourier Frequency f 

Figure 5. Power spectral density (FFT) of the error term of Mangoldt function A(n). 

Hardy found that the RFT of the modified Mangoldt function b(n) = A(n)(j)(n) jn equals 
ll(q) /(j)(q). It is thus interesting to look at the summatory function 

B(t) = £ A(n)*(n)/n = t{\ +£ B (t)). (29) 

n=\ 

The error term in Fig. 6 is found to follow approximately the power law 

s B {t)^r 2a oo) 

with a = (v5— l)/2 = 1/(1 + 1/(1 + 1/[1 + ...))), as shown in Fig. 7. This spectrum shows 
a possible connection between a and /i(g) and thus a possible relationship between the theory 
of diophantine approximations for quadratic irrational numbers such as a and prime number 
theory. The RFT of £ B (t) looks similar to the one n(q)/(j)(q) of the new Mangoldt function 
b{n). 
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2000 4000 6000 8000 10000 12000 14000 16000 
time t 



Figure 6. Error term in the new Mangoldt function b{n). 




1 10 100 1000 

Fourier Frequency f 



Figure 7. Power spectral density (FFT) of the error term in new Mangoldt function b(n) in 
comparison to the power law 1 /f 2a , with a = (y5 — l)/2, the golden mean. 




5 10 15 20 25 

resonance q 



Figure 8. Ramanujan-Fourier transform (RFT) of the error term (upper curve) of new 
Mangoldt function b(n) in comparison to the function jj.(q) /<p(q)(lower curve). 

4. LOW FREQUENCY NOISE FROM EXPERIMENTAL DATA 

Our final goal in using the Ramanujan-Fourier transform is to discover known arithmetical 
rules behind experimental sequences. 

4.1. Low frequency noise from galactic nuclei 



We give an example taken from astronomy. The observation of variability in astronomical 
systems may lead to valuable information on the physical nature of the observed system. In 
particular Seyfert galaxies are a subset of galaxies which exhibit evidence for highly energetic 
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phenomena in their nuclei: they are called active galactic nuclei or AGN. They are thought to 
be powered by accretion onto massive black holes at their centers. X-rays are created mainly 
in high temperature, high density regimes, and since matter is fairly transparent to high energy 
X-rays, monitoring X-ray emission from AGNs provides a view into the core and may be used 
to understand the accretion process there. 

Here we used a sample of data taken from the EXOSAT archive by M. Koenig and 
available at http://astro.uni-tuebingen.de/groups/time/ (sample mkn766-85.dat. outZRM) (see 
Fig. 9). The power spectral density exhibits a 1/f low frequency noise as well as white noise 
as shown in Fig. 10. The corresponding RFT analysis shown in Fig. 1 1 shows a well defined 
signature reminiscent of the RFT signature of Mangoldt function, that is n(q)/(j)(q). That 
may be an indication that many resonance processes occur between the black hole and the 
matter to be accreted, a process which may be described from prime number theory. 




500 1000 1500 2000 2500 3000 3500 4000 
time t 



Figure 9. X ray variability from an active galactic nucleus (AGN). 



0.1 




1 10 100 1000 

Fourier Frequency f 

Figure 10. Power spectral density (FFT) of X ray variability from an AGN. 



4.2. Low frequency noise close to phase locking 

Our second example is taken from the study of radio frequency oscillators close to phase 
locking. We recently demonstrated a relation between phase locking, \/f frequency noise 
and prime numbers [6]. According to that approach the coupling coefficient between the 
oscillators could be described from a Mangoldt function, leading to de-synchronization effects 
and 1/f frequency noise. The RFT should be able to support that conjecture. Figs. 12, 13 and 
14 show the beat note close to phase locking of 5 MHz oscillators, the 1/ / noise calculated 
from the FFT and the corresponding RFT. 
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Figure 11. Ramanujan-Fourier transform (RFT) of X ray variability of an AGN. 
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Figure 12. Beat frequency between two radio frequency oscillators close to phase locking. 
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Figure 13. FFT of the beat frequency for two oscillators close to phase locking. 
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Figure 14. RFT of the beat frequency for two oscillators close to phase locking. 



5. DISCUSSION 



There are at least two great challenges in the study of Ramanujan sums. One can be interested 
in the extraction of arithmetical features from experimental files, with the aim to develop a 
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relevant theory of their randomness. We have in mind the signals exhibiting 1/ f noise, since 
this type of noise still carries much mystery, in electronics as well as in other fields, from 
physics to biology and society. 

The RFT signature of 1/ f noise in a phase locked loop studied in Sect. 4.2 still doesn't 
keep one's promise, since we were unable to relate it to a known arithmetical function. 
Further work is required. In contrast, high energy astrophysics seems to be a relevant field for 
Ramanujan sums based signal processing. They may help to derive plausible theories of the 
strong variability observed close to Seyfert or other massive galaxies. See the reference [12] 
for another application of arithmetic to the black-hole remote sensing problem. 

The other challenge behind Ramanujan sums relates to prime number theory. We just 
focused our interest to the relation between 1/ / noise in communication circuits and the still 
unproved Riemann hypothesis[13]. The mean value of the modified Mangoldt function b(n), 
introduced in (29), links Riemann zeros to the 1/ f 2a noise and to the Mobius function. This 
should follow from generic properties of the modular group 5L(2,Z), the group of 2 by 2 
matrices of determinant 1 with integer coefficients [14], and to the statistical physics of Farey 
spin chains[15]. See also the link to the theory of Cantorian fractal space time[16]. 
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